Calling  in  sick:  impacts  of  fever  on  intra-urban  human  mobility 
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Abstract 

Pathogens  inflict  a  wide  variety  of  disease  manifestations  on  their  hosts,  yet  the  impacts  of 
disease  on  the  behavior  of  infected  hosts  are  rarely  studied  empirically  and  are  seldom  accounted 
for  in  mathematical  models  of  transmission  dynamics.  We  explored  the  potential  impacts  of  one 
of  the  most  common  disease  manifestations,  fever,  on  a  key  detenninant  of  pathogen 
transmission,  host  mobility,  in  residents  of  the  Amazonian  city  of  Iquitos,  Peru.  We  did  so  by 
comparing  two  groups  of  febrile  individuals  (dengue-positive  and  dengue-negative)  with  an 
afebrile  control  group.  A  retrospective,  semi-structured  interview  allowed  us  to  quantify  multiple 
aspects  of  mobility  during  the  two- week  period  preceding  each  interview.  We  fitted  nested 
models  of  each  aspect  of  mobility  to  data  from  interviews  and  compared  models  using  likelihood 
ratio  tests  to  detennine  whether  there  were  statistically  distinguishable  differences  in  mobility 
attributable  to  fever  or  its  etiology.  Compared  to  afebrile  individuals,  febrile  study  participants 
spent  more  time  at  home,  visited  fewer  locations,  and,  in  some  cases,  visited  locations  closer  to 
home  and  spent  less  time  at  certain  types  of  locations.  These  multifaceted  impacts  are  consistent 
with  the  possibility  that  disease-mediated  changes  in  host  mobility  generate  dynamic  and 
complex  changes  in  host  contact  network  structure. 

Background 

It  is  well  known  that  pathogens  are  capable  of  effecting  a  wide  range  of  disease  manifestations  in 
their  hosts.  Disease  manifestations  can  significantly  alter  behavior  in  human  and  non-human 
animal  hosts  [1],  leading  to  important  consequences  for  the  ecology  of  both  host  [2]  and 
pathogen  [3]  populations.  They  can  also  mediate  behavioral  changes  in  susceptible,  uninfected 


hosts  seeking  to  avoid  infection  and  illness  [4].  Such  behavioral  responses  by  susceptible  hosts 
have  been  recently  shown  to  exert  important  dynamic  feedbacks  on  the  course  of  epidemics  [4- 
6], 

One  effect  of  many  disease  manifestations  on  host  behavior  is  modification  of  host 
mobility;  e.g.,  changes  in  home  range,  distance  and  speed  of  movement,  or  locations  visited. 
Although  pathogens  such  as  Toxoplasma  gondii  have  been  shown  to  greatly  increase  host 
mobility  in  rodents  [7],  it  is  likely  that  disease  manifestations  more  typically  lead  to  an  overall 
reduction  in  mobility — anyone  who  has  ever  stayed  home  from  work  or  school  due  to  illness  is  a 
case  in  point.  Such  effects  have  potentially  broad  significance  throughout  epidemiology, 
ecology,  and  evolution.  Hypothesized  consequences  of  disease-mediated  reductions  in  mobility 
range  from  shifts  in  the  seasonality  of  avian  influenza  outbreaks  in  waterfowl  [8]  to  the  evolution 
of  virulence  attenuation  in  directly  transmitted  pathogens  as  compared  with  vector-bome 
pathogens,  such  as  Plasmodium  spp.  [9],  More  generally,  how  disease  impacts  host  mobility  and 
contact  has  been  proposed  as  a  key  determinant  of  how  within-host  infection  dynamics  affect 
population-level  transmission  dynamics  [10]. 

Despite  the  expected  ubiquity  of  disease  impacts  on  mobility  and  their  documentation  in 
animal  systems  [11],  they  have  seldom  been  measured  empirically  in  humans.  One  of  the  few 
studies  of  this  phenomenon  in  humans  was  based  on  self-completed  written  questionnaires 
administered  by  mail  to  179  Britons  diagnosed  with  swine  flu  [12].  That  study,  which  focused  on 
the  measurement  of  contacts  suitable  for  influenza  transmission,  showed  that  illness  reduced 
contacts  by  approximately  50%.  A  follow-up  study  [13]  estimated  that  reductions  in  contact  due 
to  the  impacts  of  disease  on  behavior  caused  a  71%  reduction  in  the  pathogen’s  basic 
reproductive  number,  R0,  compared  to  what  it  would  have  been  had  sick  people  engaged  in  as 
many  contacts  as  healthy  people.  Thus,  it  is  clear  that  the  pervasive  effects  of  disease 
manifestations  on  host  behavior  could  have  important  epidemiological  consequences,  yet  these 
effects  are  almost  completely  absent  from  transmission  models  of  infectious  diseases  of  humans 
(although  see  [14]). 

Quantifying  the  impacts  of  illness  on  direct  contacts  [12,13]  was  a  significant  advance  for 
understanding  the  epidemiology  of  influenza  and  other  directly  transmitted  pathogens. 
Nonetheless,  there  remain  a  number  of  gaps  in  knowledge  about  the  impacts  of  illness  on  contact 
that  merit  additional  study.  First,  for  indirectly  transmitted  pathogens  and  for  exposure  to 
environmentally  mediated  agents  of  disease,  a  highly  relevant  concept  is  that  of  the  activity 
space  [15,16],  which  is  a  description  of  a  person’s  time  allocation  across  the  set  of  all  locations 
that  she  or  he  frequents  [17],  similar  to  the  concept  of  a  utilization  distribution  in  the  animal 
movement  literature  [18].  Establishing  that  illness  reduces  direct  contacts  does  not  address  how 
illness  might  affect  time  allocation  at  key  locations  for  exposure  to  environmental  contaminants 
or  biting  by  insect  vectors.  Second,  a  gap  that  could  not  possibly  be  filled  by  the  aforementioned 
study  [12]  alone  is  that  it  pertained  to  a  single  geographic  and  cultural  context,  leaving  open  the 
possibility  of  different  behavioral  consequences  of  illness  under  different  circumstances.  People 
in  resource-poor,  tropical  environments,  for  example,  are  known  to  have  comparatively  complex 
patterns  of  intra-urban  mobility,  which  has  important  consequences  for  transmission  dynamics 
[19].  Third,  the  nature  and  severity  of  disease  manifestations  vary  considerably  across  different 
diseases  and  could,  therefore,  impact  behavior,  mobility,  and  contact  in  different  ways. 

To  begin  to  address  these  gaps  in  the  understanding  of  disease  impacts  on  mobility,  we 
used  retrospective,  semi-structured  interviews  of  950  participants  in  a  longitudinal  study  of 
dengue  epidemiology  in  the  Amazonian  city  of  Iquitos,  Peru.  Because  these  data  came  from  a 
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longitudinal  study  of  dengue,  a  disease  caused  by  viruses  transmitted  by  day-biting  mosquitoes, 
interviews  were  designed  to  characterize  individuals’  activity  spaces  rather  than  their  contacts 
with  other  people  [20,21].  Symptoms  associated  with  dengue  are  classically  described  as  acute 
fever,  headache,  musculoskeletal  pain,  and  rash  [22,23],  For  analytical  tractability,  however,  we 
restricted  our  analysis  to  a  binary  classification  of  study  participants  as  either  febrile  or  afebrile 
based  on  a  threshold  body  temperature  of  38  °C.  Nonetheless,  there  is  extensive  variation  in  the 
severity  of  fever  and  other  symptoms  of  dengue,  which  allowed  disease  to  impact  multiple 
features  of  the  measured  activity  spaces  of  our  study  participants.  Additionally,  while  some 
febrile  participants  had  laboratory  confirmed  dengue  virus  infections  (DENV+),  others  did  not 
(DENV-)  and  were  instead  likely  infected  with  some  other  pathogen,  such  as  influenza  or  a 
bacterial  agent.  This  allowed  us  to  explore  the  possibility  that  impacts  of  fever  on  mobility 
depend  on  the  fever’s  etiology. 

Methods 

Overview 

Our  primary  goal  was  to  identify  aspects  of  human  mobility  for  which  there  are  clear  and 
unambiguous  effects  of  fever.  To  do  so,  we  used  interview  data  collected  from  residents  of 
Iquitos,  Peru  (salient  characteristics  of  which  are  described  in  the  Supplemental  Methods). 
Interview  data  were  stratified  by  demographic  status,  fever  status,  and  fever  etiology  to  perform 
separate  comparisons  among  these  strata  for  each  aspect  of  mobility  that  we  considered.  We 
performed  comparisons  under  a  hypothesis  testing  framework;  e.g.,  is  a  given  aspect  of  mobility 
perceptibly  different  between  febrile  and  afebrile  people?  We  prioritized  the  identification  of 
effects  over  detailed  quantification  of  effect  sizes  or  probabilistic  statements  about  what  level  of 
detail  best  explains  the  data.  Although  effect  sizes  or  probabilistic  statements  could  be  insightful, 
we  restricted  our  analysis  to  a  series  of  hypothesis  tests  because  of  strengths  and  limitations  of 
our  interview  data  and  because  of  the  need  to  establish  basic  knowledge  about  the  impacts  of 
fever  on  human  mobility. 

Mobility  data  collection  and  processing 

We  collected  data  on  intra-urban  mobility  by  characterizing  study  participants’  activity  spaces 
with  retrospective,  semi-structured  movement  interviews  (SSI).  The  SSI  was  designed  in 
accordance  with  focus  group  discussions  with  local  residents  [20],  validated  against  comparable 
data  collected  using  Global  Positioning  System  data-loggers  [21],  and  used  to  parameterize 
models  of  intra-urban  human  mobility  [29].  A  detailed  description  of  the  SSI  is  available  in  the 
Supplemental  Methods. 

At  its  conclusion,  each  interview  yielded  a  table  describing  the  activity  space  of  a  study 
participant  in  the  two-week  period  prior  to  the  interview.  Each  row  of  this  table  corresponded  to 
a  location  visited  by  the  study  participant,  and  respective  columns  of  this  table  contained  a 
unique  identifier  for  each  location,  its  latitude,  its  longitude,  its  designation  as  belonging  to  one 
of  the  eight  aforementioned  land  use  types,  the  frequency  per  day  of  visits  to  the  location  by  the 
study  participant,  and  the  average  duration  of  each  visit  to  the  location  by  the  study  participant. 
These  tables  were  then  used  to  derive  tables  that  aggregated  specific  features  of  the  activity 
spaces  of  all  study  participants,  including:  (1)  a  table  containing  the  total  number  of  locations 
visited  by  each  study  participant;  (2)  eight  tables  containing  the  number  of  locations  of  each  type 
visited  by  each  study  participant;  (3)  eight  tables  containing  the  geographic  coordinates  of  each 
location  of  a  given  type  that  was  visited  by  a  given  study  participant,  as  well  as,  the  geographic 
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coordinates  of  that  study  participant’s  home;  (4)  nine  tables  containing  the  frequency  per  day  of 
visits  to  each  study  participant’s  home  and  to  locations  of  each  type  visited  by  each  study 
participant;  and  (5)  nine  tables  containing  the  average  duration  of  visits  to  each  study 
participant’s  home  and  to  locations  of  each  type  visited  by  each  study  participant. 

Study  participant  enrollment  and  classification 

Participants  in  our  study  were  enrolled  as  part  of  a  longitudinal,  cluster-based  investigation  of 
dengue  epidemiology.  In  total,  that  study  involved  2,444  participants  aged  5  years  or  older  living 
in  two  intensively  studied  neighborhoods,  Maynas  and  Tupac  Amaru  (Fig.  SI).  This  longitudinal 
cohort  study  has  been  described  in  full  detail  by  Stoddard  et  al.  [26].  Salient  details  of  that  study 
include  that  febrile  participants  were  identified  by  members  of  the  field  team  who  actively 
monitored  for  fever  (defined  as  body  temperature  >38  °C)  in  each  house,  which  was  visited  at 
least  three  times  per  week  during  the  study  period.  Blood  was  drawn  upon  detection  of  fever 
(acute  sample)  and  approximately  15  days  later  (convalescent  sample),  and,  for  participants  who 
were  willing,  an  SSI  was  performed  within  24  hrs  of  detection  of  fever  in  92%  of  cases.  In 
addition,  residents  of  other  houses  visited  by  febrile  participants  were  recruited  into  the  study 
and  asked  to  participate  in  an  SSI  upon  initiation,  rather  than  requiring  detection  of  fever.  We 
designated  study  participants  who  were  interviewed  without  a  detectable  fever  as  “afebrile”  and 
those  who  participated  in  an  SSI  at  the  time  at  which  fever  was  detected  as  “febrile.” 

In  addition  to  being  febrile  or  afebrile,  we  further  classified  study  participants  according 
to  two  dimensions  (Table  1).  First,  febrile  participants  were  classified  as  either  having  an  active 
dengue  infection  at  the  time  of  the  SSI  (DENV+)  or  having  a  fever  for  presumably  some  other 
reason  (DENV-).  This  determination  was  made  by  testing  acute  and  convalescent  blood  samples 
for  dengue-specific  IgM  antibodies  by  antibody-capture  ELISA  and  by  testing  acute  samples  for 
dengue  viruses  by  real-time  polymerase  chain  reaction  [31],  Second,  we  categorized  study 
participants  according  to  four  demographic  classifications  that  we  presume  account  for  a  large 
portion  of  inter-individual  variation  in  mobility.  Classifications  included:  school-age  children 
aged  5-17  (whose  mobility  is  likely  dominated  by  school,  home,  and  neighbors’  houses),  college 
students  aged  >18  (whose  mobility  is  likely  dominated  by  attending  college  but  may  have  more 
variability  than  children),  homemakers  and  unemployed  adults  (who  presumably  spend  more 
time  at  home  than  other  adults),  and  adults  aged  >18  who  work  outside  the  home  (who  likely 
spend  a  great  deal  of  time  outside  the  home  at  diverse  locations). 

The  high  degree  of  infonnality  in  the  Iquitos  labor  market  [19]  precluded  further 
refinement  of  the  working  adult  population,  because  a  large  proportion  of  the  population  has 
more  than  one  job,  leading  to  over  150  self-reported  occupations  among  our  study  participants. 
Although  impacts  of  fever  on  mobility  were  our  primary  interest,  accounting  for  other  sources  of 
variation  in  mobility  was  important,  because  it  reduced  the  potential  for  confounding;  e.g.,  there 
were  relatively  more  school-age  children  among  febrile  participants  than  there  were  among 
afebrile  participants  (Table  1).  Had  we  not  disaggregated  study  participants  according  to  these 
demographic  categories,  differences  in  movement  attributable  solely  to  demographic  differences 
could  have  been  mistakenly  attributed  to  differences  associated  with  fever  or  its  etiology. 

Models  and  statistical  analyses 

For  each  aspect  of  mobility  encoded  in  the  aforementioned  tables,  we  used  models  with  the  same 
level  of  detail  and  parametric  forms  used  in  the  analysis  by  Perkins  et  al.  [29],  with  some  minor 
modifications.  For  the  number  of  locations  visited,  we  assumed  a  negative  binomial  distribution 
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for  the  total  number  of  locations  of  all  types,  but  a  Poisson  distribution  for  the  number  of 
locations  of  each  type.  We  modeled  the  probability  of  choosing  a  location  of  a  given  distance 
from  home  d  as  the  product  of  the  number  of  locations  at  distance  d  weighted  by  an  exponential 
function  exp(-p  d)  with  a  unique  parameter  p  for  each  location  type.  For  the  frequency  and 
duration  of  periods  of  time  at  home  and  at  locations  of  each  type,  we  assumed  separate  bivariate 
normal  distributions  on  a  log-log  scale  with  non-zero  correlations,  and  we  assumed  that  distance 
from  home  had  no  effect  on  the  frequency  or  duration  of  visits. 

To  test  hypotheses  about  differences  in  mobility  attributable  to  fever,  its  etiology,  and 
demographic  category,  we  fitted  models  for  each  aspect  of  mobility  with  study  participants 
agglomerated  in  different  ways.  For  example,  in  the  simplest  model,  a  given  aspect  of  mobility 
for  all  study  participants  was  assumed  to  be  described  by  a  single  set  of  parameters.  By  contrast, 
in  the  most  complex  model,  a  given  aspect  of  mobility  was  described  by  distinct  parameter  sets 
for  each  of  twelve  groups:  four  demographic  categories  by  three  fever  statuses  (i.e.,  febrile  and 
DENV+,  febrile  and  DENV-,  afebrile).  With  respect  to  fever  status,  we  considered  models 
denoted  F3  (mobility  of  each  group  explained  by  distinct  parameters),  F2  (mobility  of 
individuals  with  no  fever  distinct  from  those  with  any  fever),  and  FI  (mobility  indistinguishable 
with  respect  to  fever  status).  With  respect  to  demographic  category,  we  considered  models 
denoted  C4  (mobility  of  individuals  in  each  group  explained  by  distinct  parameters)  and  C 1 
(mobility  of  individuals  indistinguishable  with  respect  to  demographic  category).  The  letter  in 
each  of  these  codes  denotes  the  factor  (i.e.,  F  =  febrile  status,  C  =  demographic  category), 
whereas  the  numeral  denotes  the  number  of  distinct  subgroups  considered  in  the  corresponding 
model;  e.g.,  2  =  each  individual  fell  into  one  of  two  possible  categories.  Combining  models  of 
fever  status  and  demographic  category  yielded  the  following  set  of  models  describing  both: 
F3,C4;  F2,C4;  F1,C4;  F3,C1;  and  F2,C1. 

To  compare  these  models,  and  thereby  to  test  hypotheses  about  the  influence  of  fever 
status  and  demographic  category  on  each  aspect  of  mobility,  we  perfonned  a  series  of  pairwise 
likelihood  ratio  tests  for  each  aspect  of  mobility.  Likelihood  ratio  tests  are  applicable  to  nested 
models  and  quantify  the  probability  that  the  observed  deviance  between  two  models  is  due  to 
chance  [32].  For  example,  F1,C4  is  nested  within  F2,C4,  because  the  former  is  a  special  case  of 
the  latter,  but  with  some  of  the  parameters  constrained  to  be  equal.  For  each  aspect  of  mobility, 
we  fitted  all  candidate  models  by  likelihood  maximization  and  performed  a  likelihood  ratio  test 
for  each  pair  of  nested  models  in  the  R  language  [33].  Given  twelve  pairs  of  nested  models  for 
each  aspect  of  mobility,  we  addressed  each  hypothesis  by  considering  the  results  of  those 
statistical  tests  on  the  whole  rather  than  by  parsing  each  statistical  test  as  corresponding  to  a 
meaningful  and  distinct  hypothesis.  For  example,  if  a  difference  in  mobility  attributable  to  fever 
was  detected  when  all  participants  were  pooled,  but  not  when  they  were  disaggregated  by 
demographic  category  (e.g.,  F2,C1  >  FI, Cl  but  F2,C4  =  F1,C4),  then  our  overall  conclusion  was 
that  there  was  no  difference  in  mobility  attributable  to  fever  because  the  difference  that  was 
detected  could  be  attributable  solely  to  demographic  differences.  Furthermore,  because  we 
performed  twelve  separate  tests  for  each  aspect  of  mobility,  we  regarded  support  for  significance 
between  two  models  as  strongest  if  the  /7-value  of  a  test  was  less  than  a=0. 05/12=4. 17xl0"3, 
which  was  obtained  by  applying  a  Bonferoni  correction  [34].  Although  our  analysis  could  have 
been  performed  under  alternative  frameworks,  such  as  infonnation  theoretic  or  Bayesian 
approaches,  our  choice  of  a  series  of  nested  likelihood  ratio  tests  was  appropriate  given  our 
hypothesis  testing  goals  and  the  nested  relationships  among  our  models  [35]. 
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Results 

Fever  had  consistent  and  discernable  effects  on  all  aspects  of  mobility  that  we  examined 
(summarized  in  Table  2).  For  the  total  number  of  locations  that  people  visited,  febrile  study 
participants  visited  30%  fewer  locations  on  average  than  afebrile  participants  (Fig.  1),  with  all 
pairwise  tests  of  models  with  F3  or  F2  against  models  with  FI  having  p  <  10"13  (Table  SI).  This 
difference  could  be  due  to  febrile  participants  visiting  fewer  locations  of  many  types,  but  given 
the  available  data  and  the  statistical  approach  that  we  used,  commercial  locations  were  the  only 
location  type  for  which  there  were  consistent  reductions  across  all  demographic  categories  (Fig. 
S5).  All  pairwise  tests  of  models  with  F3  or  F2  against  FI  had  p  <  1CT10  (Table  S3).  This  effect 
was  weaker  for  school-age  children,  but  for  other  demographic  categories,  febrile  participants 
visited  at  least  one  fewer  commercial  location  on  average  (Fig.  S5). 

There  were  also  differences  in  the  tendencies  for  febrile  and  afebrile  participants  to  visit 
locations  of  varying  distances  from  home.  After  removing  the  effect  of  the  distribution  of 
locations  within  a  given  distance  of  a  person’s  home,  parameters  describing  the  strength  of  one’s 
preference  for  visiting  locations  closer  to  home  differed  significantly  between  febrile  and  afebrile 
participants  and  between  those  who  were  DENV+  and  DENV-,  with  all  pairwise  tests  of  models 
with  F3  against  F2  or  FI  having p  <  10"4  (Table  S10).  We  could  not  identify  systematic 
differences  (e.g.,  we  could  not  conclude  that  febrile  participants  consistently  exhibited  a  stronger 
preference  for  visiting  locations  closer  to  home),  but  instead  detected  wholesale  differences 
between  multiple  sets  of  eight  location-type-specific  parameters  (Fig.  S12).  Although  not  tested 
statistically  on  an  individual  basis,  effects  that  appeared  particularly  notable  and  that  seemed 
logical  included  that  febrile,  DENV+  adults  who  do  not  work  outside  the  home  tended  to  restrict 
their  movements  to  residential  and  commercial  locations  closer  to  home  than  their  afebrile 
counterparts  did,  and  school-age  children  and  college  students  visited  recreational  locations 
closer  to  home  when  febrile  (Figs.  S12  &  S13). 

Febrile  study  participants  generally  spent  more  time  at  home  than  afebrile  participants, 
although  the  extent  of  this  effect  varied  across  demographic  categories  (Fig.  2).  Febrile  school- 
age  children  spent  26.6%  more  time  at  home  than  their  afebrile  counterparts  when  DENV+  and 
21.6%  more  when  DENV-.  Among  adults,  those  who  do  not  work  outside  the  home  exhibited  the 
greatest  increase  in  time  spent  at  home  due  to  fever  (DENV+:  3.8%  more  than  afebrile,  DENV-: 
12.8%  more  than  afebrile),  whereas  differences  in  time  spent  at  home  between  febrile  and 
afebrile  college  students  and  adults  working  outside  the  home  were  less  consistent  (Fig.  2).  A 
pairwise  test  of  F2,C1  and  FI, Cl  models  was  statistically  significant  (p  =  0.003),  as  were  tests  of 
F2,C4  against  F1,C4  (p  <  10"12)  and  tests  comparing  models  with  F3  against  those  with  F2  (p  < 
10'3)  (Table  SI  1).  The  only  location  type  besides  home  for  which  there  was  a  significant 
difference  in  time  spent  per  location  was  the  educational  location  type.  Febrile  school-age 
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children  spent  less  time  at  each  educational  location  they  visited  than  afebrile  children  did, 
primarily  due  to  a  difference  in  the  frequency  of  visits  rather  than  the  duration  of  visits  (Fig. 

SI  7).  Pairwise  tests  of  models  with  differences  in  both  fever  status  and  demographic  category 
were  statistically  significant  with p  <  0.01  (Table  S15). 

Beyond  differences  attributable  to  fever  on  the  whole,  there  were  no  clear  differences  in 
any  aspect  of  mobility  between  febrile  study  participants  with  differing  etiologies.  The  majority 
of  pairwise  comparisons  of  models  with  F3  against  those  with  F2  were  not  significant  (Tables 
SI -SI 9),  and  those  that  were  significant  tended  to  support  differences  that  were  relatively  weak 
or  equivocal  compared  to  differences  attributable  to  fever  on  the  whole  (Figs.  1,  2,  S12).  For 
example,  although  there  were  statistically  significant  differences  in  the  tendency  to  visit 
locations  closer  to  home  among  afebrile,  DENV+  febrile,  and  DENV-  febrile  study  participants 
( p  <  1CT3,  Table  S10),  there  were  no  clear,  systematic  differences  across  all  location  types  (Figs. 
S12&S13). 

Irrespective  of  fever  status,  there  were  substantial  differences  in  mobility  among 
participants  from  different  demographic  categories.  School-age  children  tended  to  visit  fewer 
unique  locations  overall  (Fig.  1),  spend  more  time  at  home  (Fig.  2),  visit  locations  closer  to  home 
(Figs.  S12  &  S13),  visit  more  residential  locations  (Fig.  S4),  and  spend  more  time  at  each 
residential  (Fig.  S14)  and  educational  location  that  they  visited  (Fig.  S17)  than  participants  from 
the  other  demographic  categories.  College  students  also  spent  more  time  at  educational  locations 
that  they  visited  (Fig.  SI 7),  but  were  otherwise  similar  to  other  adults.  Adults  who  did  not  work 
outside  the  home  were  distinguished  by  spending  more  time  at  home  than  participants  from  other 
categories  (Fig.  2). 

There  also  appeared  to  be  a  select  few  interactions  between  the  effects  of  fever  and 
demographic  category  on  mobility.  School-age  children  tended  to  visit  approximately  the  same 
number  of  commercial  locations  (2. 5 -2. 7)  regardless  of  their  febrile  status,  whereas  febrile 
participants  from  other  demographic  categories  consistently  visited  on  average  one  fewer 
commercial  location  when  febrile  (Fig.  S5,  Table  S3).  Differences  between  participants’  overall 
preferences  for  visiting  locations  closer  to  home  also  depended  on  demographic  category,  as 
stated  earlier  (Fig.  S12,  Table  S10).  Time  spent  at  home  by  college  students  did  not  appear  to 
depend  on  fever  status,  whereas  febrile  participants  from  other  demographic  categories  spent 
significantly  more  time  at  home  than  their  afebrile  counterparts  (Fig.  2). 

Discussion 

We  used  retrospective,  semi-structured  interviews  to  assess  a  very  basic,  but  surprisingly 
neglected,  question  about  host-pathogen  interaction  in  humans:  what  is  the  impact  of  fever  on  an 
infected  person’s  mobility?  Our  results  revealed  differences  in  all  aspects  of  mobility  between 
study  participants  who  were  febrile  during  the  time  they  were  interviewed  and  those  who  were 
afebrile.  Reductions  in  the  number  of  locations  visited  and  time  spent  outside  the  home  were 
among  the  most  unambiguous  differences  between  febrile  and  afebrile  participants,  regardless  of 
their  demographic  category.  These  and  other  effects,  such  as  visiting  locations  closer  to  home 
and  changes  in  the  types  of  locations  visited,  are  logically  consistent  with  one  another  and  with 
previous  reports  from  urban  settings  in  developed  countries  [12,13].  These  results  are  also 
consistent  with  the  concept  of  “sickness  behaviors”  in  the  animal  behavior  literature,  where 
impaired  mobility  and  other  effects  are  better  studied  and  are  thought  to  have  important  adaptive 
consequences  [1 1,36].  Together,  our  results  suggest  that  fever  significantly  impairs  human 
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mobility  and,  by  extension,  could  modulate  the  potential  for  infectious  people  with  clinically 
apparent  infections  to  engage  in  contacts  and  contribute  to  onward  transmission. 

The  effects  that  we  observed  pertained  to  a  two-week  period  prior  to  presentation  of 
fever,  which  exceeds  the  duration  of  fever  for  people  infected  with  dengue,  influenza,  and  many 
other  pathogens  [37,38].  Thus,  it  is  most  likely  that  participants  in  our  study  experienced  fever 
for  only  a  few  days  prior  to  their  interviews.  If  so,  that  would  imply  that  impacts  of  fever  on 
mobility  during  the  time  when  study  participants  were  actually  febrile  may  have  been  stronger 
than  the  effects  that  we  were  able  to  measure.  For  example,  rather  than  spending  26.6%  more 
time  at  home  (Fig.  2),  this  would  suggest  that  DENV+  febrile  children  may  have  spent  nearly  all 
of  their  time  at  home  on  days  when  they  experienced  fever.  Similarly,  some  effects  could  have 
gone  undetected  altogether.  For  example,  it  is  noteworthy  that  we  detected  no  increase  in  visits 
to  healthcare  locations  among  febrile  individuals.  This  is  not  altogether  surprising,  however, 
given  the  common  occurrence  of  relatively  mild  disease  associated  with  dengue  infection  [39]. 
This  result  could  also  be  an  artifact  of  our  study  design,  as  some  study  participants  may  have 
waited  for  a  visit  from  a  project  physician  rather  than  seeking  healthcare  at  a  local  clinic  or 
hospital.  For  these  and  similar  reasons,  our  quantitative  results  should  be  interpreted  as 
conservative  underestimates  of  the  impacts  of  fever  on  mobility. 

Although  our  results  are  of  limited  value  quantitatively,  they  provide  an  important 
starting  point  for  determining  the  extent  to  which  the  impacts  of  disease  manifestations  on 
mobility  are  perceptible  with  interviews,  one  of  several  common  instruments  for  measuring 
detailed,  individual-level  mobility  [40,41],  One  limitation  of  the  retrospective  SSI  that  we  used  is 
that  it  relied  on  study  participant  recall,  which  is  subject  to  recall  bias.  Even  so,  in  a  previous 
study  in  Iquitos  [21],  we  compared  data  obtained  from  the  SSI  to  comparable  data  obtained  from 
wearable  GPS  data-loggers  and  found  that,  overall,  the  SSI  was  a  superior  means  of  measuring 
activity  spaces  due  to  behavioral  limitations  associated  with  GPS  methods;  e.g.,  forgetting  to 
wear  the  GPS  data-logger.  With  regard  to  this  study,  it  is  possible  that  recall  bias  may  have 
impaired  our  ability  to  detect  the  full  extent  of  the  true  effect  of  fever  on  mobility,  in  which  case 
our  estimates  of  the  impact  of  fever  on  mobility  are  conservative.  This  could  be  so  if  individuals 
who  have  visited  more  locations  (i.e.,  afebrile  individuals)  are  less  likely  to  recall  all  the 
locations  they  have  visited  compared  to  those  who  have  visited  fewer  locations;  i.e.,  febrile 
individuals.  Addressing  this  and  other  open  questions — such  as  the  impact  of  other  disease 
manifestations  on  mobility  or  how  suppression  of  fever  with  antipyretics  affects  mobility  during 
key  times  in  the  course  of  a  person’s  infectiousness — was  beyond  the  scope  of  our  study,  but 
constitutes  an  important  direction  for  future  studies  designed  with  similar  objectives  in  mind. 

Future  studies  should  examine  behavioral  modification  due  to  fever  across  the  period  of 
time  that  a  person  is  infectious.  For  dengue,  infectiousness  often  lasts  longer  than  fever.  In 
experimental  feeding  studies,  infection  of  mosquitoes  was  documented  up  to  two  days  before 
[42]  and  two  days  after  [43]  the  period  of  fever,  which  typically  lasts  between  two  and  seven 
days  [37].  In  a  more  recent  study  [44],  at  a  given  viremia  level,  dengue  virus  infected  people 
with  no  detectable  symptoms  or  whose  symptoms  had  not  yet  begun  were  more  infectious  to 
mosquitoes  than  people  who  were  bitten  when  they  were  symptomatic.  Fully  addressing  how 
fever  and  other  disease  symptoms  affect  an  individual’s  contribution  to  transmission  will  require 
combining  data  on  human  mobility  with  data  on  infectiousness  at  multiple  time  points  over  the 
course  of  infection  [10,45]. 

Pathogen  natural  history,  particularly  mode  of  transmission,  will  also  affect  how  impacts 
of  fever  on  mobility  influence  impacts  of  fever  on  contact  and  opportunities  for  pathogen 
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transmission.  For  dengue  and  other  vector-bome  diseases,  the  impact  of  spending  more  time  at 
home  and  visiting  fewer  locations  on  vector-human  contact  will  depend  on  the  densities  of 
competent  vectors  at  those  locations  [15].  If  there  are  more  vectors  at  a  person’s  home  than  at  the 
locations  where  they  spend  less  time  due  to  fever,  then  they  could  effectively  be  engaging  in 
more  mosquito  contacts  than  if  fever  had  not  reduced  their  mobility.  Conversely,  if  there  are 
fewer  vectors  at  home  than  elsewhere  in  a  person’s  activity  space,  the  converse  would  be  true. 
Densities  of  the  DENV  mosquito  vector  Aedes  aegypti  in  Iquitos  are  sufficiently  heterogeneous 
in  time  and  space  that  both  of  these  possibilities  probably  happen  [43].  For  directly  transmitted 
pathogens,  spending  more  time  at  home  and  visiting  fewer  locations  should  result  in  fewer 
opportunities  for  transmission,  on  average  [13].  Other  effects  of  fever  could  counteract  this 
tendency,  including  increased  visitation  by  members  of  a  sick  person’s  family  or  social  network 
[6],  Furthermore,  because  the  impacts  of  fever  on  an  infected  person’s  mobility  are  inherently 
dynamic  and,  as  we  have  shown,  multifaceted,  they  are  likely  to  generate  dynamic  and  complex 
changes  in  contact  network  structure  over  time.  For  both  directly  and  indirectly  transmitted 
pathogens,  mathematical  modeling  will  play  an  essential  role  in  future  studies  seeking  to 
elucidate  how  disease-mediated  changes  in  the  mobility  of  infected  hosts  affect  the  dynamic 
nature  of  contact  networks  and  pathogen  transmission. 

Accurately  quantifying  a  realistic  balance  between  the  various  effects  of  fever  on  contact 
and  the  relationship  between  fever  and  infectiousness  is  important  for  applications  of  population- 
level  models  to  forecast  future  dynamics,  infer  information  about  key  parameters  from  historical 
dynamics,  and  assess  intervention  impacts  [5,46].  In  addition  to  disease-mediated  changes  in 
contact  rates  within  a  population,  it  is  also  likely  that  fever  and  other  disease  manifestations 
could  modify  mobility  and  contact  between  humans  and  mosquitoes.  Spatial  models  used  for 
forecasting  the  spread  of  pandemics  [47]  and  for  estimating  sources  and  sinks  of  transmission 
[48]  routinely  utilize  movement  data  from  flight  records,  call  data  records,  and  other  sources,  or 
a  continuum  of  models  derived  from  physical  principles  [49],  such  as  gravity  [50]  and  radiation 
[51]  models.  These  approaches  rely  on  data  or  assumptions  about  healthy  individuals  and  do  not 
consider  how  disease  manifestations  might  impact  mobility  and,  consequently,  spatial 
transmission  dynamics.  Based  on  effects  that  we  detected — such  as  a  tendency  for  febrile  people 
to  visit  fewer  locations  and  to  spend  more  time  at  home — spatial  models  ignoring  these  effects 
could  overestimate  the  potential  for  a  disease  to  become  pandemic  or  overemphasize  the 
importance  of  mobility  in  identifying  sources  and  sinks  of  transmission.  Better  understanding  the 
impacts  of  disease  manifestations  on  mobility  during  times  when  hosts  are  infectious  merits 
increased  attention  in  infectious  disease  epidemiology. 
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Table  1.  Numbers  of  study  participants  who  provided  information  about  time  spent  at  home  or 
time  spent  at  locations  other  than  home,  stratified  by  fever  status  and  demographic  category. 


School-age 

children 

College  student 

Homemaker  / 
Unemployed 
adult 

Working  adult 

Time  at  home 

Afebrile 

34 

14 

36 

44 

Febrile,  DENV+ 

91 

15 

24 

41 

Febrile,  DENY- 

214 

53 

94 

116 

Time  elsewhere 

Afebrile 

43 

25 

52 

53 

Febrile,  DENV+ 

143 

25 

32 

59 

Febrile,  DENY- 

241 

66 

121 

128 
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Table  2.  Summary  of  differences  between  febrile  and  afebrile  study  participants  with  respect  to 
different  aspects  of  mobility. 


Aspect  of  mobility 

Impact  of  fever 

Tables  and  Figures 

Total  number  of 
locations  visited 

Febrile  study  participants  from  all  demographic 
categories  consistently  visited  fewer  locations. 

Table  SI,  Figure  1 

Types  of  locations 
visited 

Febrile  study  participants  from  all  demographic 
categories  consistently  visited  fewer  commercial 
locations.  Consistent  differences  between  the 
numbers  of  locations  of  all  other  types  that  febrile 
and  afebrile  participants  visited  were  not  as 
perceptible. 

Table  S3,  Figure  S5; 
Tables  S2  &  S4-S9, 
Figures  S4&S6-S11 

Distance  from 

On  the  whole,  there  were  significant  differences 

Table  S10,  Figures 

home  of  locations 

between  the  distances  from  home  to  locations 

S12&S13 

visited 

visited  by  febrile  and  afebrile  study  participants, 
particularly  when  accounting  for  demographic 
category.  Notable  differences  included  that 
homemakers  with  dengue  fever  visited  residential 
and  commercial  locations  closer  to  home,  and 
febrile  children  and  college  students  visited 
recreational  locations  closer  to  home. 

Time  spent  at 
home 

Febrile  study  participants  from  all  demographic 
categories  consistently  spent  more  time  at  home 
overall,  which  resulted  from  coming  and  going 
from  home  less  frequently  and  from  staying  there 
longer  each  time  they  were  home. 

Table  S12,  Figure  2 

Time  spent  per 
location  at 

The  most  perceptible  difference  in  time  spent  per 
location  at  locations  other  than  home  was  that 

Table  S15,  Figure 

S17;  Tables  S12-S14 

locations  other 

febrile  school-aged  children  spent  less  time  at 

&  S16-S19,  Figures 

than  home 

educational  locations.  This  difference  appeared  to 
be  driven  more  by  a  lower  frequency  of  visits  than 
by  a  shorter  duration. 

S14-S16  &  S18-S21 
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Figure  1.  Number  of  locations  visited  by  study  participants,  stratified  by  participant  fever  status 
(rows)  and  demographic  category  (columns).  Gray  bars  show  the  empirical  distribution,  the 
black  line  shows  the  empirical  mean  (numerical  value  printed  in  each  panel),  and  the  balls  and 
stems  show  the  fitted  negative  binomial  distribution  with  maximum-likelihood  estimates  of 
parameter  values  (size  and  prob  as  in  dnbinom  function  in  R,  values  indicated  in  each  panel). 
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Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  2.  Joint  probability  distributions  of  the  frequency  and  mean  duration  of  periods  of  time 
spent  at  home.  Colors  represent  continuous  probabilities  that  range  low  (red),  medium  (orange), 
and  high  (bright  yellow  and  white).  Parameter  values  of  each  distribution  were  fitted  by 
likelihood  maximization  to  data  from  individuals  from  each  of  four  groups  (children,  students, 
homemakers,  others)  using  interviews  conducted  when  those  individuals  either  had  no  fever, 
dengue  fever,  or  fever  due  to  some  other  cause.  Sample  means  of  duration,  frequency,  and 
proportion  of  total  time  spent  at  home  for  each  group  are  noted  in  each  panel. 


School-age  children 


College  students 


Homemakers  /  Unemployed  adults 


Working  adults 


Frequency  of  visits  (times  per  day) 
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Supplemental  Tables  and  Figures 

Table  SI.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
locations  visited. 

Table  S2.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
residential  locations  visited. 

Table  S3.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
commercial  locations  visited. 

Table  S4.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
recreational  locations  visited. 

Table  S5.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
educational  locations  visited. 

Table  S6.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
religious  locations  visited. 

Table  S7.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
healthcare  locations  visited. 

Table  S8.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
institutional  locations  visited. 

Table  S9.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
locations  of  other  types  visited. 

Table  S10.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  the 
distance  from  home  of  locations  visited  by  study  participants. 

Table  Sll.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  home. 

Table  S12.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  residential  locations  other  than  home. 

Table  S13.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  commercial  locations. 

Table  S14.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  recreational  locations. 

Table  S15.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  educational  locations. 
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Table  S16.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  healthcare  locations. 

Table  S17.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  religious  locations. 

Table  S18.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  institutional  locations. 

Table  S19.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  locations  of  some  type  other  than  those  considered  previously. 

Figure  SI.  Map  of  Iquitos  showing  the  two  neighborhoods  in  which  study 
participants  resided. 

Figure  S2.  Dates  on  which  semi-structured  interviews  took  place. 

Figure  S3.  Age  and  sex  distribution  of  study  participants. 

Figure  S4.  Number  of  residential  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category. 

Figure  S5.  Number  of  commercial  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category. 

Figure  S6.  Number  of  recreational  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category. 

Figure  S7.  Number  of  educational  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category. 

Figure  S8.  Number  of  religious  locations  visited  by  study  participants,  stratified  by 
participant  fever  status  and  demographic  category. 

Figure  S9.  Number  of  healthcare  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category. 

Figure  S10.  Number  of  institutional  locations  visited  by  study  participants, 
stratified  by  participant  fever  status  and  demographic  category. 

Figure  Sll.  Number  of  locations  of  other  types  visited  by  study  participants, 
stratified  by  participant  fever  status  and  demographic  category. 

Figure  S12.  Maximum-likelihood  estimates  (MLE)  of  the  parameter  for  the  function 
that  weights  distance  from  home  on  one’s  probability  of  visiting  a  location  of  a  given 
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distance  from  home,  stratified  by  participant  fever  status  and  demographic 
category. 


Figure  S13.  Mean  distance  from  home  of  locations  visited  by  study  participants, 
stratified  by  participant  fever  status  and  demographic  category. 

Figure  S14.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  residential  locations  other  than  home. 

Figure  SI  5.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  commercial  locations. 

Figure  S16.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  recreational  locations. 

Figure  S17.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  educational  locations. 

Figure  S18.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  healthcare  locations. 

Figure  S19.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  religious  locations. 

Figure  S20.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  institutional  locations. 

Figure  S21.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  locations  of  some  type  other  than  those  considered  previously. 
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Table  SI.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
locations  visited.  The  p-value  associated  with  each  test  is  the  probability  of 
obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value.  Parameter 
values  for  each  model  were  fitted  by  likelihood  maximization  of  a  negative  binomial 
distribution  of  the  number  of  locations  that  each  study  participant  visited. 
Maximum-likelihood  estimates  of  parameter  values  for  participants  with  each 
combination  of  fever  status  and  demographic  category  are  shown  in  Fig.  1. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xl  =  25.1 

Xi6  =  144.6 

Xis  =  61.2 

xlo  =  82.1 

X22  =  145.0 

p  =  0.001 

p<  io-15 

p<  10-5 

p  <  IO-8 

p<  10-15 

F2,C4 

— 

xl  =  119.5 
p<  io-15 

— 

X12  =  57.0 

p  <  io-7 

X14  =  H9.8 

p<  10-15 

F1,C4 

— 

— 

— 

— 

xl  =  0.34 
p  =  0.99 

F3,C1 

— 

— 

— 

xl  =  20.9 

P  <  IO'4 

Xl  =  83.8 

p<  10-15 

F2,C1 

— 

— 

— 

— 

Xl  =  62.8 

p<  10-13 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S2.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
residential  locations  visited.  The  p-value  associated  with  each  test  is  the  probability 
of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value. 
Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  residential  locations  that  each  study 
participant  visited  other  than  her  or  his  own  home.  Maximum-likelihood  estimates 
of  parameter  values  for  participants  with  each  combination  of  fever  status  and 
demographic  category  are  shown  in  Fig.  S4. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

Xa  =  12.0 

xl  =  19.8 

Xl  =  21.3 

xlo  =  24.2 

xli  =  25.9 

p  =  0.017 

p  =  0.011 

p  =  0.012 

p  =  0.007 

p  =  0.007 

F2,C4 

— 

Xa  =  7.8 

p  =  0.10 

— 

xl  =  12.2 

p  =  0.057 

xl  =  13.9 
p  =  0.054 

F1,C4 

— 

— 

— 

— 

Xl  =  6-0 

p  =  0.11 

F3,C1 

— 

— 

— 

Xl  =  2.9 
p  =  0.087 

xl  =  4.6 

p  =  0.10 

F2,C1 

— 

— 

— 

— 

xl  =  1.6 

p  =  0.19 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S3.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
commercial  locations  visited.  The  p-value  associated  with  each  test  is  the  probability 
of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value. 
Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  commercial  locations  that  each  study 
participant  visited.  Maximum-likelihood  estimates  of  parameter  values  for 
participants  with  each  combination  of  fever  status  and  demographic  category  are 
shown  in  Fig.  S5. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

N 

II 

xl  =  67.1 

Xl  =  23.5 

Xio  =  28.1 

Xli  =  86.9 

p  =  0.10 

p<  io-10 

p  =  0.005 

p  =  0.002 

P<  10-!3 

F2,C4 

— 

xl  =  59.4 
p<  io-11 

— 

Xl  =  20.4 

p  =  0.002 

Xi  =  79.2 

P<  10-!3 

F1,C4 

— 

— 

— 

— 

CO 

■  00 

2s 

II  V 

F3,C1 

— 

— 

— 

xl  =  4.6 
p  =  0.032 

xl  =  63.4 
p<  lO-w 

F2,C1 

— 

— 

— 

— 

Zl2  =  58.8 
p<  lO-^ 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S4.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
recreational  locations  visited.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  recreational  locations  that  each  study 
participant  visited.  Maximum-likelihood  estimates  of  parameter  values  for 
participants  with  each  combination  of  fever  status  and  demographic  category  are 
shown  in  Fig.  S6. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xl  =  5-8 

Xl  =  13-2 

Xl  =  7-4 

Xw  =  9.7 

*i2i  =  14.1 

p  =  0.22 

p  =  0.10 

p  =  0.60 

p  =  0.47 

p  =  0.23 

F2,C4 

— 

Xa  =  7.5 

p  =  0.11 

— 

Xl  =  3-9 
p  =  0.69 

Xi  =  8.3 
p  =  0.30 

F1,C4 

— 

— 

— 

— 

Xl  =  0.9 
p  =  0.83 

F3,C1 

— 

— 

— 

Xi2  =  2.3 
p  =  0.13 

Xl  =  6-7 
p  =  0.035 

F2,C1 

— 

— 

— 

— 

Zi2  =  4.4 
p  =  0.036 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S5.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
educational  locations  visited.  The  p-value  associated  with  each  test  is  the  probability 
of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value. 
Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  educational  locations  that  each  study 
participant  visited.  Maximum-likelihood  estimates  of  parameter  values  for 
participants  with  each  combination  of  fever  status  and  demographic  category  are 
shown  in  Fig.  S7. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

X4  =  2.2 

Xl  =  3-1 

Xl  =  3-5 

xlo  =  4.4 

*i2i  =  4.9 

p  =  0.69 

p  =  0.93 

p  =  0.94 

p  =  0.93 

p  =  0.93 

F2,C4 

— 

Xl  =  0.8 
p  =  0.93 

— 

Xl  =  2-2 
p  =  0.91 

xl  =  2.7 
p  =  0.91 

F1,C4 

— 

— 

— 

— 

Xl  =  1.9 

p  =  0.60 

F3,C1 

— 

— 

— 

Xl  =  0.9 
p  =  0.34 

*22  =  1.4 
p  =  0.50 

F2,C1 

— 

— 

— 

— 

Xl  =  0.5 
p  =  0.48 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S6.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
religious  locations  visited.  The  p-value  associated  with  each  test  is  the  probability  of 
obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value.  Parameter 
values  for  each  model  were  fitted  by  likelihood  maximization  of  a  Poisson 
distribution  of  the  number  of  religious  locations  that  each  study  participant  visited. 
Maximum-likelihood  estimates  of  parameter  values  for  participants  with  each 
combination  of  fever  status  and  demographic  category  are  shown  in  Fig.  S8. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xl  =  0-47 

Xl  =  0.99 

Xl  =  0.79 

xio  =  i.o 

xh  =  i.o 

p  =  0.98 

p  =  0.99 

p  =  0.99 

p  =  0.99 

p  =  0.99 

F2,C4 

— 

Xl  =  0-5 
p  =  0.97 

— 

Xl  =  0-55 
p  =  0.99 

Xi  =  0.6 
p  =  0.99 

F1,C4 

— 

— 

— 

— 

xi  =  0.04 
p  =  0.99 

F3,C1 

— 

— 

— 

Xi2  =  0.2 
p  =  0.62 

xi  =  0.02 

p  =  0.89 

F2,C1 

— 

— 

— 

— 

Zi2  <  10-5 
p  =  0.99 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S7.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
healthcare  locations  visited.  The  p-value  associated  with  each  test  is  the  probability 
of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value. 
Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  healthcare  locations  that  each  study 
participant  visited.  Maximum-likelihood  estimates  of  parameter  values  for 
participants  with  each  combination  of  fever  status  and  demographic  category  are 
shown  in  Fig.  S9. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xt  =  0.24 

Xl  =  1-2 

Xl  =  1-7 

Xto  =  1-9 

xti  =  2.6 

p  =  0.99 

p  =  0.99 

p  =  0.99 

p  =  0.99 

p  =  0.99 

F2,C4 

— 

xl  =  1-0 

p  =  0.91 

— 

Xl  =  1-7 
p  =  0.95 

Xi  =  2.4 
p  =  0.94 

F1,C4 

— 

— 

— 

— 

Xl  =  1-4 
p  =  0.71 

F3,C1 

— 

— 

— 

Xl  =  0-2 
p  =  0.64 

X4  =  0.9 
p  =  0.64 

F2,C1 

— 

— 

— 

— 

=  0.7 
p  =  0.41 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S8.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
institutional  locations  visited.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  institutional  locations  that  each  study 
participant  visited.  Maximum-likelihood  estimates  of  parameter  values  for 
participants  with  each  combination  of  fever  status  and  demographic  category  are 
shown  in  Fig.  S10. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xl  =  1-9 

Xl  =  4-5 

Xl  =  6.7 

xlo  =  7.5 

*i2i  =  7.9 

p  =  0.75 

p  =  0.81 

p  =  0.67 

p  =  0.68 

p  =  0.72 

F2,C4 

— 

X4  =  2.5 
p  =  0.64 

— 

Xl  =  5-6 
p  =  0.47 

xl  =  6.0 

p  =  0.54 

F1,C4 

— 

— 

— 

— 

Xl  =  3-5 
p  =  0.32 

F3,C1 

— 

— 

— 

Xl  =  0.8 
p  =  0.36 

xl  =  1-3 
p  =  0.53 

F2,C1 

— 

— 

— 

— 

Xl  =  0.4 
p  =  0.51 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S9.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  number  of 
locations  of  other  types  visited.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
Poisson  distribution  of  the  number  of  locations  of  other  types  that  each  study 
participant  visited.  Maximum-likelihood  estimates  of  parameter  values  for 
participants  with  each  combination  of  fever  status  and  demographic  category  are 
shown  in  Fig.  Sll. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xl  =  0-2 

Xl  =  5-8 

Xl  =  3.7 

X10  =  3.7 

xh  =  10.6 

p  =  0.99 

p  =  0.67 

p  =  0.93 

p  =  0.96 

p  =  0.47 

F2,C4 

— 

x\  =  5.7 

p  =  0.22 

— 

Xl  =  3-6 
p  =  0.73 

xl  =  10.5 

p  =  0.16 

F1,C4 

— 

— 

— 

— 

Xl  =  4-8 
p  =  0.19 

F3,C1 

— 

— 

— 

Ji2  <  IQ"3 
p  =  0.98 

*22  =  6.9 
p  =  0.03 

F2,C1 

— 

— 

— 

— 

Zi2  =  6.9 

p  =  0.01 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S10.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  the 
distance  from  home  of  locations  visited  by  study  participants.  The  p-value 
associated  with  each  test  is  the  probability  of  obtaining  a  deviance  between  Ha  and 
Ho  greater  than  the  observed  value.  Parameter  values  for  each  model  were  fitted  by 
likelihood  maximization  of  exponential  weighting  functions  for  the  effect  of  distance 
from  home  on  the  probability  of  choosing  a  location  of  a  given  type  and  distance 
from  home.  Maximum-likelihood  estimates  of  parameter  values  for  participants 
with  each  combination  of  fever  status  and  demographic  category  are  shown  in  Fig. 
S12. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

Xh  =  68.9 

X264  =  H6.3 

X72  =  766.4 

Xio  =  801.3 

*88  =  836.7 

P<  lO-3 

P<  lO'4 

p<  1015 

p<  1015 

p<  1015 

F2,C4 

— 

X32  =  47.4 
p  =  0.039 

— 

X43  =  732.5 

p<  10-15 

Xs6  =  767.9 

p<  1015 

F1,C4 

— 

— 

— 

— 

Zl4  =  720.4 

p<  1015 

F3,C1 

— 

— 

— 

Xl  =  35.0 

P  <  lO'4 

X16  =  70.4 

p  <  10-8 

F2,C1 

— 

— 

— 

— 

Xl  =  35.4 

P<  lO'4 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  Sll.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  home.  The  p-value  associated  with  each  test  is  the  probability  of  obtaining  a 
deviance  between  Ha  and  H0  greater  than  the  observed  value.  Parameter  values  for 
each  model  were  fitted  by  likelihood  maximization  of  a  bivariate  lognormal 
distribution  of  frequency  and  mean  duration  of  periods  of  time  spent  at  home.  The 
mean  time  spent  at  home,  as  well  as  the  mean  frequency  and  duration  of  bouts  at 
home,  for  participants  with  each  combination  of  fever  status  and  demographic 
category  is  shown  in  Fig.  4. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

X20  =  55.0 

Xio  =  158.1 

Xls  =  250.5 

Xio  =  274.6 

Xh  =  292.5 

p<  10-4 

p<  10-15 

p<  1015 

p<  1015 

p<  1015 

F2,C4 

— 

xio  =  103.1 

p<10-12 

— 

Xio  =  219.6 

p<  1015 

Xls  =  237.6 

p<  1015 

F1,C4 

— 

— 

— 

— 

xh  =  134.4 

p<  1015 

F3,C1 

— 

— 

— 

Xl  =  24.1 
p<  10-3 

Xio  =  42.0 
p  <  lO-5 

F2,C1 

— 

— 

— 

— 

Xl  =  17.9 
p  =  0.003 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S12.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  residential  locations  other  than  home.  The  p-value  associated  with  each  test 
is  the  probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the 
observed  value.  Parameter  values  for  each  model  were  fitted  by  likelihood 
maximization  of  a  bivariate  lognormal  distribution  of  frequency  and  mean  duration 
of  periods  of  time  spent  at  residential  locations  other  than  home.  The  mean  time 
spent  at  each  residential  location,  as  well  as  the  mean  frequency  and  duration  of 
visits  to  those  locations,  for  participants  with  each  combination  of  fever  status  and 
demographic  category  is  shown  in  Fig.  S14. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

*20  =  27.8 

*lo  =  47.9 

*45  =  90-7 

Xlo  =  104.8 

*55  =  112.2 

p  =  0.11 

p  =  0.18 

P  <  10-4 

p<  10-5 

p<  10-5 

F2,C4 

— 

*20  =  20.1 
p  =  0.45 

— 

Xlo  =  77.0 
p  <  10-5 

*35  =  84.4 
p  <  10-5 

F1,C4 

— 

— 

— 

— 

*15  =  64.3 

p<  io-7 

F3,C1 

— 

— 

— 

*5  =  14.1 
p  =  0.015 

*i2o  =  21.6 
p  =  0.017 

F2,C1 

— 

— 

— 

— 

*5  =  7.5 
p  =  0.19 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S13.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  commercial  locations.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
bivariate  lognormal  distribution  of  frequency  and  mean  duration  of  periods  of  time 
spent  at  commercial  locations.  The  mean  time  spent  at  each  residential  location,  as 
well  as  the  mean  frequency  and  duration  of  visits  to  those  locations,  for  participants 
with  each  combination  of  fever  status  and  demographic  category  is  shown  in  Fig. 
S15. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

lio  =  14.4 

Xlo  =  77.7 

xls  =  69-3 

Xlo  =  69-8 

Xls  =  119-4 

p  =  0.81 

P<  lO-3 

p  =  0.012 

p  =  0.033 

p<  lO'5 

F2,C4 

— 

X20  =  63.3 
p  <  lO  5 

— 

xlo  =  55.4 
p  =  0.003 

Xls  =  105.0 
p  <  lO-8 

F1,C4 

— 

— 

— 

— 

Xh  =  41.7 

P<  lO'3 

F3,C1 

— 

— 

— 

xl  =  0.59 
p  =  0.99 

X10  =  50.1 
p  <  lO  6 

F2,C1 

— 

— 

— 

— 

xl  =  49.5 

p  <  10'8 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S14.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  recreational  locations.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
bivariate  lognormal  distribution  of  frequency  and  mean  duration  of  periods  of  time 
spent  at  recreational  locations.  The  mean  time  spent  at  each  residential  location,  as 
well  as  the  mean  frequency  and  duration  of  visits  to  those  locations,  for  participants 
with  each  combination  of  fever  status  and  demographic  category  is  shown  in  Fig. 
S16. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

Xlo  =  9-2 

X40  =  62.4 

xls  =  55.7 

Xlo  =  60.3 

Xls  =  H7.6 

p  =  0.98 

p  =  0.013 

p  =  0.13 

p  =  0.15 

p<  10-5 

F2,C4 

— 

xlo  =  53.2 

p  <  10-4 

— 

Xlo  =  51.1 
p  =  0.009 

Xls  =  108.4 
p  <  lO-8 

F1,C4 

— 

— 

— 

— 

Xis  =  55.2 
p<  lO'5 

F3,C1 

— 

— 

— 

xl  =  4.6 
p  =  0.46 

xlo  =  61.9 

p  <  10-8 

F2,C1 

— 

— 

— 

— 

xl  =  57.2 
p  <  10-10 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S15.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  educational  locations.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
bivariate  lognormal  distribution  of  frequency  and  mean  duration  of  periods  of  time 
spent  at  educational  locations.  The  mean  time  spent  at  each  residential  location,  as 
well  as  the  mean  frequency  and  duration  of  visits  to  those  locations,  for  participants 
with  each  combination  of  fever  status  and  demographic  category  is  shown  in  Fig. 
S17. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

X20  =  30.0 

xio  =  66.7 

xls  =  277.7 

Xio  =  293.4 

xis  =  327.2 

p  =  0.07 

p  =  0.005 

p  <  10-15 

p<  1015 

p<  1015 

F2,C4 

— 

xio  =  36.8 

p  =  0.012 

— 

Xio  =  263.5 

p<  1015 

Xh  =  297.3 

p<  1015 

F1,C4 

— 

— 

— 

— 

Xi5  =  260.5 

p<  1015 

F3,C1 

— 

— 

— 

Xl  =  15.7 

p  =  0.008 

xio  =  49.5 
p  <  lO-6 

F2,C1 

— 

— 

— 

— 

Xl  =  33.8 
p<  10-5 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S16.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  healthcare  locations.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
bivariate  lognormal  distribution  of  frequency  and  mean  duration  of  periods  of  time 
spent  at  healthcare  locations.  The  mean  time  spent  at  each  residential  location,  as 
well  as  the  mean  frequency  and  duration  of  visits  to  those  locations,  for  participants 
with  each  combination  of  fever  status  and  demographic  category  is  shown  in  Fig. 
S18. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

X20  =  15.3 

xlo  =  51.9 

xls  =  110-8 

Xlo  =  H5.4 

Xls  =  154.4 

p  =  0.76 

p  =  0.10 

p  <  lO-6 

p  <  lO-6 

p  <  1010 

F2,C4 

— 

Xlo  =  36.6 
p  =  0.013 

— 

Xlo  =  100.1 

p  <  IO-8 

Xls  =  139.2 

p<1013 

F1,C4 

— 

— 

— 

— 

X15  =  102.6 
p  <  IO-14 

F3,C1 

— 

— 

— 

xl  =  4.59 
p  =  0.47 

Xio  =  43.6 

p  <  io-5 

F2,C1 

— 

— 

— 

— 

Xl  =  39.0 
p  <  IO-6 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S17.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  religious  locations.  The  p-value  associated  with  each  test  is  the  probability 
of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed  value. 
Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
bivariate  lognormal  distribution  of  frequency  and  mean  duration  of  periods  of  time 
spent  at  religious  locations.  The  mean  time  spent  at  each  residential  location,  as  well 
as  the  mean  frequency  and  duration  of  visits  to  those  locations,  for  participants  with 
each  combination  of  fever  status  and  demographic  category  is  shown  in  Fig.  S19. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

X20  =  20.2 

xlo  =  33.6 

Xls  =  42-6 

Xlo  =  50.5 

xls  =  53.3 

p  =  0.44 

p  =  0.75 

p  =  0.58 

p  =  0.45 

p  =  0.54 

F2,C4 

— 

Xlo  =  13.3 

p  =  0.86 

— 

X30  =  30.3 
p  =  0.45 

Xls  =  33.0 
p  =  0.56 

F1,C4 

— 

— 

— 

— 

Xis  =  19.7 

p  =  0.18 

F3,C1 

— 

— 

— 

Xl  =  7.96 

p  =  0.16 

Xio  =  10.7 
p  =  0.38 

F2,C1 

— 

— 

— 

— 

Xl  =  2-7 
p  =  0.74 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 


20 


Table  S18.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  institutional  locations.  The  p-value  associated  with  each  test  is  the 
probability  of  obtaining  a  deviance  between  Ha  and  H0  greater  than  the  observed 
value.  Parameter  values  for  each  model  were  fitted  by  likelihood  maximization  of  a 
bivariate  lognormal  distribution  of  frequency  and  mean  duration  of  periods  of  time 
spent  at  institutional  locations.  The  mean  time  spent  at  each  residential  location,  as 
well  as  the  mean  frequency  and  duration  of  visits  to  those  locations,  for  participants 
with  each  combination  of  fever  status  and  demographic  category  is  shown  in  Fig. 
S20. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

*20  =  27.8 

xlo  =  68.9 

*45  =  75-7 

Xlo  =  80.0 

xis  =  100.6 

p  =  0.12 

p  =  0.003 

p  =  0.003 

p  =  0.005 

p<  lO-3 

F2,C4 

— 

Xlo  =  41.2 
p  =  0.004 

— 

Xlo  =  52.2 
p  =  0.007 

*35  =  72-8 

P  <  lO'3 

F1,C4 

— 

— 

— 

— 

*15  =  31.6 
p  =  0.007 

F3,C1 

— 

— 

— 

xl  =  4.3 
p  =  0.51 

*i2o  =  24.9 

p  =  0.006 

F2,C1 

— 

— 

— 

— 

*5  =  20.6 

P  <  lO'3 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Table  S19.  Results  from  likelihood  ratio  tests  for  pairs  of  nested  models  of  time 
spent  at  locations  of  some  type  other  than  those  considered  previously.  The  p-value 
associated  with  each  test  is  the  probability  of  obtaining  a  deviance  between  Ha  and 
Ho  greater  than  the  observed  value.  Parameter  values  for  each  model  were  fitted  by 
likelihood  maximization  of  a  bivariate  lognormal  distribution  of  frequency  and 
mean  duration  of  periods  of  time  spent  at  locations.  The  mean  time  spent  at  each 
residential  location,  as  well  as  the  mean  frequency  and  duration  of  visits  to  those 
locations,  for  participants  with  each  combination  of  fever  status  and  demographic 
category  is  shown  in  Fig.  S21. 


Ha 

Ho 

F2,C4 

F1,C4 

F3,C1 

F2,C1 

FI, Cl 

F3,C4 

xlo  =  20.6 

Xlo  =  46.7 

xls  =  66.5 

Xlo  =  78.9 

Xls  =  97.4 

p  =  0.42 

p  =  0.22 

p  =  0.02 

p  =  0.006 

p<  lO-3 

F2,C4 

— 

X\  o  =  26.1 

p  =  0.16 

— 

xlo  =  58.3 

p  =  0.001 

Xls  =  76.8 

p  <  10-4 

F1,C4 

— 

— 

— 

— 

Xis  =  50.6 
p<  lO'5 

F3,C1 

— 

— 

— 

Xl  =  12.5 
p  =  0.029 

Xio  =  30.9 
p  <  lO'3 

F2,C1 

— 

— 

— 

— 

Xl  =  18.4 
p  =  0.003 

F3  =  3  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile  &  DENV+,  febrile  &  DENV-) 

F2  =  2  different  classes  with  respect  to  fever  (i.e.,  afebrile,  febrile) 

FI  =  no  distinction  between  movement  of  afebrile  and  febrile 

C4  =  4  different  classes  of  people  (i.e.,  school-age  children,  college  student,  homemaker  /  unemployed  adult,  working  adult) 
Cl  =  no  distinction  among  people  based  on  demographic  category 
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Latitude 


Figure  SI.  Map  of  Iquitos  showing  the  two  neighborhoods  in  which  study 
participants  resided.  The  Maynas  neighborhood  is  shown  in  blue,  Tupac  Amaru  is 
shown  in  red,  and  the  location  of  Iquitos  is  indicated  by  a  yellow  star. 
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Figure  S2.  Dates  on  which  semi-structured  interviews  took  place,  stratified  by  the 
fever  status  of  study  participants  across  the  rows  and  binned  by  month.  The  left 
(right)  column  includes  participants  who  participated  in  the  semi-structured 
interview  pertaining  to  time  spent  at  home  (locations  other  than  home). 


Time  at  home  Time  elsewhere 


2008  1  2009  1  2010  2011  >  2012  >  2013  2008  >  2009  i  2010  2011  i  2012  >  2013 


2008  2009  2010  2011  2012  2013  2008  2009  2010  2011  2012  2013 


Date 


24 


Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  S3.  Age  and  sex  distributions  of  study  participants  stratified  by  fever  status 
along  the  columns.  The  top  (bottom)  row  includes  participants  who  participated  in 
the  semi-structured  interview  pertaining  to  time  spent  at  home  (locations  other 
than  home). 


Afebrile  Febrile,  DENV+  Febrile,  DENV- 


cr 

CD 


Age 


25 


Away  Home 


Figure  S4.  Number  of  residential  locations  other  than  home  visited  by  study 
participants,  stratified  by  participant  fever  status  and  demographic  category.  Gray 
bars  show  the  empirical  distribution,  the  black  line  shows  the  empirical  mean,  and 
the  balls  and  stems  show  the  fitted  Poisson  distribution  with  maximum-likelihood 
estimates  of  parameter  values.  The  mean  number  of  locations  of  this  type  that  were 
visited  is  displayed  in  the  bottom  right  of  each  panel. 
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Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  S5.  Number  of  commercial  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category.  Gray  bars  show  the  empirical 
distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and  stems  show 
the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of  parameter 
values.  The  mean  number  of  locations  of  this  type  that  were  visited  is  displayed  in 
the  bottom  right  of  each  panel. 
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Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  S6.  Number  of  recreational  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category.  Gray  bars  show  the  empirical 
distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and  stems  show 
the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of  parameter 
values.  The  mean  number  of  locations  of  this  type  that  were  visited  is  displayed  in 
the  bottom  right  of  each  panel. 
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Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  S7.  Number  of  educational  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category.  Gray  bars  show  the  empirical 
distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and  stems  show 
the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of  parameter 
values.  The  mean  number  of  locations  of  this  type  that  were  visited  is  displayed  in 
the  bottom  right  of  each  panel. 
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Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  S8.  Number  of  religious  locations  visited  by  study  participants,  stratified  by 
participant  fever  status  and  demographic  category.  Gray  bars  show  the  empirical 
distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and  stems  show 
the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of  parameter 
values.  The  mean  number  of  locations  of  this  type  that  were  visited  is  displayed  in 
the  bottom  right  of  each  panel. 
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Febrile,  DENV-  Febrile,  DENV+  Afebrile 


Figure  S9.  Number  of  healthcare  locations  visited  by  study  participants,  stratified 
by  participant  fever  status  and  demographic  category.  Gray  bars  show  the  empirical 
distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and  stems  show 
the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of  parameter 
values.  The  mean  number  of  locations  of  this  type  that  were  visited  is  displayed  in 
the  bottom  right  of  each  panel. 
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Figure  S10.  Number  of  institutional  locations  visited  by  study  participants, 
stratified  by  participant  fever  status  and  demographic  category.  Gray  bars  show  the 
empirical  distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and 
stems  show  the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of 
parameter  values.  The  mean  number  of  locations  of  this  type  that  were  visited  is 
displayed  in  the  bottom  right  of  each  panel. 
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Figure  Sll.  Number  of  locations  of  other  types  visited  by  study  participants, 
stratified  by  participant  fever  status  and  demographic  category.  Gray  bars  show  the 
empirical  distribution,  the  black  line  shows  the  empirical  mean,  and  the  balls  and 
stems  show  the  fitted  Poisson  distribution  with  maximum-likelihood  estimates  of 
parameter  values.  The  mean  number  of  locations  of  this  type  that  were  visited  is 
displayed  in  the  bottom  right  of  each  panel. 
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Figure  S12.  Maximum-likelihood  estimates  (MLE)  of  the  parameter  for  the  function 
that  weights  distance  from  home  on  one’s  probability  of  visiting  a  location  of  a  given 
distance  from  home,  stratified  by  participant  fever  status  and  demographic 
category.  Higher  values  of  these  parameters  are  consistent  with  a  stronger 
preference  for  visiting  locations  closer  to  home.  The  three  fever  status  categories 
are  indicated  by  color:  blue  =  afebrile;  green  =  febrile  &  DENV+;  red  =  febrile  & 
DENV-.  The  demographic  categories  are  indicated  by  letters:  C  =  school-age 
children;  S  =  college  student;  H  =  homemaker  /  unemployed  adult;  W  =  working 
adult. 
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Figure  S13.  Mean  distance  from  home  of  locations  visited  by  study  participants, 
stratified  by  participant  fever  status  and  demographic  category.  The  three  fever 
status  categories  are  indicated  by  color:  blue  =  afebrile;  green  =  febrile  &  DENV+; 
red  =  febrile  &  DENV-.  The  demographic  categories  are  indicated  by  letters:  C  = 
school-age  children;  S  =  college  student;  H  =  homemaker  /  unemployed  adult;  W  = 
working  adult. 
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Figure  S14.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  residential  locations  other  than  home.  Colors  represent  continuous 
probabilities  that  range  low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter 
values  of  each  distribution  were  fitted  by  likelihood  maximization  to  data  from 
individuals  from  each  of  four  groups  (School-age  children,  College  student, 
Homemaker  /  Unemployed  adult,  Working  adult)  using  interviews  conducted  when 
those  individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample 
means  of  duration,  frequency,  and  proportion  of  total  time  spent  at  residential 
locations  other  than  home  for  each  group  are  noted  in  each  panel. 
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Figure  SI  5.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  commercial  locations.  Colors  represent  continuous  probabilities  that 
range  low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter  values  of  each 
distribution  were  fitted  by  likelihood  maximization  to  data  from  individuals  from 
each  of  four  groups  (School-age  children,  College  student,  Homemaker  / 
Unemployed  adult,  Working  adult)  using  interviews  conducted  when  those 
individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample  means  of 
duration,  frequency,  and  proportion  of  total  time  spent  at  commercial  locations  for 
each  group  are  noted  in  each  panel. 
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Figure  S16.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  recreational  locations.  Colors  represent  continuous  probabilities  that 
range  low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter  values  of  each 
distribution  were  fitted  by  likelihood  maximization  to  data  from  individuals  from 
each  of  four  groups  (School-age  children,  College  student,  Homemaker  / 
Unemployed  adult,  Working  adult)  using  interviews  conducted  when  those 
individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample  means  of 
duration,  frequency,  and  proportion  of  total  time  spent  at  recreational  locations  for 
each  group  are  noted  in  each  panel. 
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Figure  S17.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  educational  locations.  Colors  represent  continuous  probabilities  that 
range  low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter  values  of  each 
distribution  were  fitted  by  likelihood  maximization  to  data  from  individuals  from 
each  of  four  groups  (School-age  children,  College  student,  Homemaker  / 
Unemployed  adult,  Working  adult)  using  interviews  conducted  when  those 
individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample  means  of 
duration,  frequency,  and  proportion  of  total  time  spent  at  educational  locations  for 
each  group  are  noted  in  each  panel. 
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Figure  S18.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  healthcare  locations.  Colors  represent  continuous  probabilities  that  range 
low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter  values  of  each 
distribution  were  fitted  by  likelihood  maximization  to  data  from  individuals  from 
each  of  four  groups  (School-age  children,  College  student,  Homemaker  / 
Unemployed  adult,  Working  adult)  using  interviews  conducted  when  those 
individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample  means  of 
duration,  frequency,  and  proportion  of  total  time  spent  at  healthcare  locations  for 
each  group  are  noted  in  each  panel. 
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Figure  S19.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  religious  locations.  Colors  represent  continuous  probabilities  that  range 
low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter  values  of  each 
distribution  were  fitted  by  likelihood  maximization  to  data  from  individuals  from 
each  of  four  groups  (School-age  children,  College  student,  Homemaker  / 
Unemployed  adult,  Working  adult)  using  interviews  conducted  when  those 
individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample  means  of 
duration,  frequency,  and  proportion  of  total  time  spent  at  religious  locations  for 
each  group  are  noted  in  each  panel. 
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Figure  S20.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  institutional  locations.  Colors  represent  continuous  probabilities  that 
range  low  (red),  medium  (orange),  and  high  (bright  yellow).  Parameter  values  of  each 
distribution  were  fitted  by  likelihood  maximization  to  data  from  individuals  from 
each  of  four  groups  (School-age  children,  College  student,  Homemaker  / 
Unemployed  adult,  Working  adult)  using  interviews  conducted  when  those 
individuals  were  afebrile,  febrile  and  DENV+,  or  febrile  and  DENV-.  Sample  means  of 
duration,  frequency,  and  proportion  of  total  time  spent  at  institutional  locations  for 
each  group  are  noted  in  each  panel. 
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Figure  S21.  Joint  distributions  of  the  frequency  and  mean  duration  of  periods  of 
time  spent  at  locations  of  some  type  other  than  those  considered  previously.  Colors 
represent  continuous  probabilities  that  range  low  (red),  medium  (orange),  and  high 
(bright  yellow).  Parameter  values  of  each  distribution  were  fitted  by  likelihood 
maximization  to  data  from  individuals  from  each  of  four  groups  (School-age 
children,  College  student,  Homemaker  /  Unemployed  adult,  Working  adult)  using 
interviews  conducted  when  those  individuals  were  afebrile,  febrile  and  DENV+,  or 
febrile  and  DENV-.  Sample  means  of  duration,  frequency,  and  proportion  of  total 
time  spent  at  locations  of  some  type  other  than  those  considered  previously  for  each 
group  are  noted  in  each  panel. 
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Supplemental  Methods 

Study  area 

The  study  was  performed  in  Iquitos,  the  largest  urban  center  in  the  Peruvian  Amazon.  Human 
population  density  in  Iquitos  is  high,  with  most  of  its  -390,000  inhabitants  living  within  a  10 
km2  area  at  an  average  household  size  of  6  people  per  house  [24].  The  regional  economy  is 
informal  (33%  of  the  economically  active  population  is  either  unemployed  or  informally 
employed)  and  major  industrial  development  circumscribed  to  small-scale  fishing,  oil  and 
lumber  extraction,  and  tourism.  Since  1999,  Iquitos  has  been  the  base  of  a  long-term 
collaborative  research  study  focusing  on  the  ecological  and  epidemiological  determinants  of 
dengue  virus  transmission,  which,  in  recent  years,  incorporated  a  strong  behavioral  component 
by  focusing  on  the  role  of  human  mobility  on  dengue  transmission  [19,20,25-27].  Iquitos  can  be 
accessed  by  air,  by  boat  along  surrounding  rivers,  and  by  automobile  only  from  a  few  nearby 
towns  [28],  but  it  is  mostly  a  relatively  isolated  city.  Consequently,  by  focusing  on  intra-urban 
mobility  rather  than  movements  away  from  the  city,  we  captured  a  majority  of  the  locations  in 
study  participants’  activity  spaces  [29]. 

Semi-structured  movement  interview 

The  goal  of  the  SSI  was  to  capture  positional  and  temporal  information  on  routine  human 
mobility  of  study  participants  [20].  It  consisted  of  a  series  of  questions  to  trigger  recall  about  the 
locations  that  the  interviewee  visited  during  the  two  weeks  preceding  the  interview,  as  well  as 
the  frequency  and  duration  of  visits  to  each  location.  To  allow  for  geocoding  of  locations  visited 
by  interviewees,  the  interviewee  was  first  asked  to  provide  an  address  for  each  location.  If  the 
address  was  unknown,  the  interviewee  was  asked  to  call  someone  at  that  location  to  obtain  the 


address  using  a  mobile  phone  provided  by  the  interviewer.  If  the  address  was  still  unknown,  the 
interviewee  was  asked  to  describe  the  whereabouts  of  the  location  to  the  best  of  his  or  her  ability, 
relying  on  a  combination  of  giving  directions  to  the  location  and/or  identifying  it  on  a  printed 
map  developed  using  contemporary,  high-quality  satellite  imagery  (Quickbird,  Digitalglobe, 

CO).  In  all  cases,  a  member  of  the  field  team  would  then  visit  the  location  and  geocode  it  for 
subsequent  analyses.  The  location’s  coordinates,  along  with  an  open-ended  description  of  the 
location’s  classification  (e.g.,  friend’s  house,  butcher  shop),  were  recorded  in  a  geographic 
information  system  [30],  For  our  analyses,  we  followed  the  methods  described  by  Perkins  et  al. 
[29]  and  grouped  each  location  as  belonging  to  one  of  the  following  eight  land  use  types: 
residential,  commercial,  recreational,  educational,  religious,  healthcare,  institutional,  or  other. 


